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Abstract 

In the framework of the MSSM, we examine several simplified models where only a few super¬ 
partners are light. This allows us to study WIMP-nucleus scattering in terms of a handful of MSSM 
parameters and thereby scrutinize their impact on dark matter direct-detection experiments. Fo¬ 
cusing on spin-independent WIMP-nucleon scattering, we derive simplified, analytic expressions 
for the Wilson coefficients associated with Higgs and squark exchange. We utilize these results to 
study the complementarity of constraints due to direct-detection, flavor, and collider experiments. 
We also identify parameter configurations that produce (almost) vanishing cross sections. In the 
proximity of these so-called blind spots, we find that the amount of isospin violation may be much 
larger than typically expected in the MSSM. This feature is a generic property of parameter re¬ 
gions where cross sections are suppressed, and highlights the importance of a careful analysis of 
the nucleon matrix elements and the associated hadronic uncertainties. This becomes especially 
relevant once the increased sensitivity of future direct-detection experiments corners the MSSM 
into these regions of parameter space. 
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I. INTRODUCTION 


Establishing the microscopic nature of Dark Matter (DM) is one of the central, open 
questions in cosmology and particle physics. In the context of cold nonbaryonic DM, the 
prevailing paradigm is based on weakly interacting massive particles (WIMPs), and exten¬ 
sive theoretical and experimental resources have been devoted towards identifying viable 
candidates and developing methods to detect them. One of the most studied WIMPs arises 
in the Minimal Supersymmetric Standard Model (MSSM), where an assumed i?-parity en¬ 
sures that the lightest superpartner (LSP) is a stable neutralino x composed of bino B , wino 
W, and Higgsino H eigenstates. The mass of the LSP is expected to lie in the range of tens 
to hundreds of GeV. 

In its general form, however, the MSSM contains more than 100 parameters, most of 
which are tied to the hidden sector which breaks supersymmetry (SUSY) at some scale 
M SUS y. Since these parameters are unknown a priori , it is necessary to restrict the di¬ 
mensionality of the parameter space in order to obtain a predictive framework with which 
to undertake phenomenological analyses. One way to achieve this is to adopt a specific 
mechanism that describes high-scale SUSY-breaking in terms of a small number of parame¬ 
ters. For example, the constrained MSSM (CMSSM) with minimal supergravity [IH1] only 
involves five free parameters, but faces increasing tension I5H32] with the non-observation 
of superpartners at the LHC experiments and other observables like the measured Higgs 
mass and anomalous magnetic moment of the muon. Alternatively, one can remain agnostic 
about the features of SUSY-breaking and incorporate data-driven constraints, as in e.g. the 
p(henomenological)MSSM [GUOS], where only 19 free parameters are used to capture the 
essential features of weak-scale SUSY. 

In both approaches, long computational chains involving spectrum generators, the calcu¬ 
lation of decay rates, or the DM relic abundance are typically required in order to explore 
the relevant parameter space. This strategy has been used extensively for the CMSSM [HHI23 
and pMSSM [Tbl ITS] to analyze y-nucleon scattering and impose limits from current DM 
direct-detection experiments such as SuperCDMS [19], XENONIOO [20], and LUX [2Tj . as 
well as upcoming proposals like XENON1T [22], LUX-ZEPLIN (LZ) [23], and SuperCDMS 
SNOLAB [23] . While these parameter scans allow one to gain useful information about the 
status of the theory in light of global fits, they generally hinder attempts to clearly identify 
which contributions associated with the underlying theory parameters can have the greatest 
impact on a signal of interest. An analytical understanding of the underlying parameter 
space can instead be obtained in the context of so-called simplified models, defined 1 [30] 
to be minimal theories of weak-scale SUSY where all but a handful of the superpartners 
relevant for DM phenomenology are decoupled from the spectrum. 

For spin-independent (SI) x-nucleon scattering, the choice of simplified model is guided 
by the dominant contributions to the cross section, namely, Higgs and squark exchange (333- 
[StTjj . To date, the focus has largely concerned the role of the Higgs sector, both in the 
decoupling limit where a single SM-like Higgs h is present in the spectrum j3U| [37], or in the 


1 For a definition of “simplified models” in the context of LHC searches, see 


2 



more general case [38], [39] where the heavier CP-even Higgs II is included. This focus is 
chiefly motivated by the fact that current bounds on the masses of gluinos and (degenerate) 
squarks of the first two generations are larger than about 1 TeV m m, and so their 
contribution to the SI cross section can be safely ignored. 2 

However, the decoupling of third-generation squarks—especially stops—upsets the main 
motivation behind the MSSM, namely, its ability to stabilize the electroweak scale 
v ~ 174 GeV against loop corrections in a technically natural fashion [46fl5l] . In other 
words, if naturalness is to remain a useful criterion with which to constrain the MSSM 
parameter space, then the spectrum should (minimally) include light stops H ,2 and—due 
to SU{2) L invariance—a left-handed sbottom b L [ 451 [oTH5T] . While the search for top and 
bottom squarks remains a primary focus of the LHC experiments, their impact on DM 
direct-detection limits has not been explored in detail. 

The purpose of this paper is to present an analytical scheme which allows one to succes¬ 
sively include those states which are most relevant for naturalness and DM direct detection. 
In particular, we consider a bino-like LSP and derive simplified, analytic expressions for the 
Wilson coefficients associated with SI scattering. We examine in detail the contributions 
from Higgs and third-generation squark exchange, and study the interplay of collider, flavor, 
and DM constraints. As in previous analyses of the Higgs sector m EHH39], our scheme 
allows us to identify so-called blind spots in parameter space, where the SI cross section is 
strongly suppressed by either a particular set of parameters mm, or destructive interfer¬ 
ence [38i [39 ] in the scattering amplitude. This effect was first identified numerically through 
a scan of the CMSSM parameter space jool I57j . while lower bounds on the yy-nucleon cross 
sections were first discussed in [58] for both the CMSSM and a generalized MSSM frame¬ 
work. A key feature of our analysis is that in the vicinity of blind spots, the amount of 
isospin violation may be much larger than typically expected for the MSSM |59]. 3 In order 
to account for isospin-violating effects originating from the nucleon matrix elements of the 
scalar quark currents, we use the formalism developed in [66], which provides an accurate 
determination of the hadronic uncertainties. 

The paper is organized as follows. In Sec. |TT| we establish our notation for y- nucleus 
scattering and comment on the differing treatments [59] E6[ m of the nucleon scalar matrix 
elements found in the literature. The leading MSSM contributions to the Wilson coefficients 
are then identified using a systematic expansion in u/Msusy which generates simplified 
analytic expressions for the Wilson coefficients associated with Higgs and squark exchange. 
Section HI examines four simplified models, where, driven by naturalness, we successively 
include the most relevant particles as active degrees of freedom. In each case, we discuss the 
conditions for blind spots and examine the amount of isospin violation allowed by current 
and projected limits from SI DM scattering. Our analysis shows that the absence of DM 


2 However, for non-degenerate squarks [121 S3] the constraints from FCNCs are satisfied [H], and the collider 
bounds are significantly weakened [45]. In this case, contributions from the first two generations could 
also be important for SI y-nucleon scattering. 

3 Large isospin violation has been put forward as a mechanism to reconcile contradictory DM direct- 
detection signal claims and null observations [50fl65| . 
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signals pushes the MSSM into regions of parameter space where isospin-violating effects are 
likely to become relevant. 


II. THEORETICAL PRELIMINARIES 

A. Spin-independent neutralino—nucleus cross section: scalar matrix elements 

We start by providing some definitions for the elastic scattering of the lightest neutralino 
X off a species of nucleus A f = %X, where Z and A denote the atomic and mass numbers 
respectively. Typically, the dependence of the cross section on the small momentum transfer 
is assumed to be described by nuclear form factors. At zero momentum transfer and for 
one-body currents only, the cross section for yA f —> yAf is given by 

4n 2 

a sl =-^[Zf p + (A-Z)f n } 2 . ( 1 ) 

7r 

Here, fi x = m x mj^/(m x + my) is the reduced mass of the y-A/" system, while f p and f n are 
effective (zero-momentum) SI couplings of the LSP to the proton and neutron respectively. 

For nucleons N, the y-A/” couplings fjy are defined by 

— = £ f? c ,+ fg £ C„, N — p or n, (2) 

q=u,d,s q=c,b,t 

where C q is the Wilson coefficient of the scalar operator fh q xx Qd with running quark mass 
fh q , and 

™ N f? = (N\rh q qq\N) , f$ = ^(1 - f» - f? - /f). (3) 

The coefficients can be interpreted as the fraction of the nucleon mass generated by the 
respective quark scalar current and are often referred to as nucleon scalar couplings. In the 
framework adopted in (J3]), the heavy quarks c,b,t are integrated out, so that, via the trace 
anomaly [681 f7T] of the QCD energy-momentum tensor, their scalar coefficients Jq can be 
expressed in terms of the light-quark ones [35]. As shown by Drees and Nojiri [72] . this 
procedure fails if the squarks are sufficiently light, and in Sec. [TT]we discuss the necessary 
modifications to ([2]) which account for the exact one-loop result. 

We note that ([3]) holds at leading order in a s : in the case of the charm quark, this may not 
be sufficiently accurate, so that either higher-order corrections [73H75] or a non-perturbative 
determination on the lattice could become mandatory. Similarly, corrections to the single¬ 
nucleon picture underlying ([Tj) in the form of two-nucleon currents can be systematically 
taken into account using effective field theory [M, 65J, [76j [77]. In this paper, we use ([TJ) 
and ([3]) to investigate the amount of isospin violation that can be generated within several 
simplified models, given the hadronic uncertainties of the single-nucleon coefficients for 
the light quarks u, d, s. 

Traditionally, the scalar matrix elements of the light quarks have been determined from a 
combination of chiral SU(3)l x SU(3)r perturbation theory (xPT 3 ) and phenomenological 
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input inferred from the pion-nucleon cr-term cr^N and the hadron mass spectrum [55, [59| 
EH; XE] • A central feature of this approach is that the up- and down-quark coefficients f/f d 
are reconstructed from two three-flavor quantities: the so-called strangeness content of the 
nucleon 


2{N\ss\N) = 2 f?/fh a 

(N\uu ± dd\N) fu/mu + fd/rh d ’ 


and another parameter 


(N\uu - ss\N) _ f/m„ - f^/m s 
(N\dd - ss\N) fd/rhd - f*/fh s 


that is related to isospin violation. As a result, the inherent uncertainties of \PT 3 (typically 
of order 30%) propagate to the two-flavor sector. Furthermore, z is usually extracted from 
a leading-order fit to baryon masses ra, and this compounds the problem of obtaining 
reliable uncertainty estimates. 

For the up- and down-quark coefficients fff d , these problems can be circumvented by 
using the two-flavor theory yPT 2 directly, thus avoiding the three-flavor expansion in the 
first place [66] . Starting from the yPT 2 expansion of the nucleon mass at third chiral order 
and including the effects due to strong isospin violation, one finds 


f N _ 

J U 

A/’ = 
A f p d = 


°±v( 1 - 0 


+ A /, 


N 


2 m N 

(1.0 ±0.2) x 10" 3 , 
(-2.1 ±0.4) x 10" 3 


fN 

Jd 

m: 

a f2 


C r ?rA r ( 1 ± 0 

2m N 


+ A/ ( 


N 


(-1.0 ±0.2) x 10" 3 , 
(2.0 ±0.4) x 10~ 3 , 


( 6 ) 


where the cr-term is defined as o^n = (N\m(uu + dd)\N), averaged over proton and neutron, 
rh = ( m u ± rhd)/ 2, and 

£ = md ~ mu = 0.36 ± 0.04 (7) 

m d ± m u 

is taken from [80] . 

For the present work, one particularly important aspect of the yPT 2 approach [66] is 
that isospin violation can be rigorously accounted for, including uncertainty estimates. This 
aspect can be nicely illustrated by considering the differences 


fu~fu = (1-9 ± 0.4) x 10 -3 and f p d - ff = (-4.1 ± 0.7) x KT 3 , (8) 

wherein the terms <t^{1 ± f)/2 from (J6]) cancel. 4 Using the yPT 3 approach, these 
differences are overestimated by roughly a factor of 2, as in e.g. R: 

/£-/»” = 4.3 xio- 3 , F d -f; = - 8.2 xio- 3 . (9) 


4 The chiral expansion of the nucleon mass difference m p — m n is known to have a large chiral logarithm 
at fourth order, with coefficient (6 g\ + l)/2 *=s 5 [HUM]. We have checked that including this logarithm 
in the analysis leads to changes of the A ff/ d well within the uncertainties given in ([6]). 
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Alternatively, one conld introduce further measures of isospin violation like and 

fu ~ fd (motivated by the quark-model picture of the nucleon), but these combinations 
depend on the specific value of cr n N. In the isospin-conserving limit, all up- and down- 
coefficients obtained from the chiral expansion of the nucleon mass become equal fu = fu = 
fd = fd = °Viv/2mAr, so that the relations /£ = and /" = are fulhllcd. 

Ultimately, the quantities relevant for the direct-detection cross section are the parameters 
defined in (j2j) , after multiplication by the Wilson coefficients and summing over quark flavors. 
In particular, the cross section Q may be rewritten as 


osi 


4 /4 f 2 

J V 


7 r 


A + (A-Z) 




( 10 ) 


so that the departure of f n j f p from unity emerges as a convenient measure of isospin viola¬ 
tion. In this context, care has to be taken in interpreting the limits on the WIMP-ra/cZeon 
cross section <7gj n given by experimental collaborations. Indeed, these are generally extracted 
via the relation 


0SI — °SI 


hx ) A 2 
Hn 


( 11 ) 


where pat is the reduced mass of the y-nucleon system. We stress that Ogj in (JTTJ) can be 
identified with the SI %-proton cross section only under the assumption f n ~ f p . If isospin- 
violating effects are large, it is natural to compare against the ^-nucleus cross section (JT|) 
directly, and indicates how the experimental limits are to be rescaled. 

In general, the effects of isospin violation depend on the target nucleus. For the mass 


range m x ~ 50-1000 GeV considered in most of our analysis (Sec. Ill), the strongest limits 
on SI X'-nucleon scattering are currently set by LUX [21]. In the context of isospin violation, 
this prompts us to focus on the projected reach of upcoming xenon-based experiments like 
XENON1T [22] and LZ [23]. However, this raises the question whether other experiments 
like SuperCDMS SNOLAB [23] (based on germanium) can be used to place complementary 
constraints on f n /f p . To quantify the difference between xenon-based constraints and other 
nuclei, consider the f n /f P dependence of the ratio 




_ KM 2 

^si y/4 e "4 Xe / 


Xe 

SI 


( 12 ) 


normalized such that Rj\r = 1 in the isospin-conserving limit. For SI scattering off argon 
and germanium, the result is shown in Fig. [lj where we observe a maximum difference of 
around 10% for f n /fp much larger or smaller than unity. In the simplified models considered 
in Sec. |III| the difference between f n and f p is generated entirely by SM quantities, so 


the improved limits offered by e.g. SuperCDMS SNOLAB are limited to the percent level. 
Moreover, the location of blind spots is determined by the condition f n>p ~ 0, so neither the 
blind spot nor the uncertainty on f n /f p depends significantly on the atomic/mass numbers 
of the nuclear target. 

The crucial input quantity oyjv is not yet precisely determined: in Sec. HI we show 
the dependence on this parameter explicitly in the case of generic Higgs exchange, and 
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fjfp 

FIG. 1: Relative difference Rj\f (defined in (|12[)) between xenon and other nuclear targets A f as 
a function of the isospin violation measure f n /fp■ The dashed line in red shows the comparison 
against N = germanium, while the solid blue line corresponds to the N = argon comparison. 


later fix its central value to cr^N = 50 MeV for illustrative purposes. The need for a precise 
determination of ovjv has triggered many ongoing efforts, including lattice-QCD calculations 
at (nearly) physical values of the pion mass; see [67J I83T185] for a compilation of recent results 
and improved phenomenological analyses. The challenge in the phenomenological approach, 
i.e. extracting oy.y from nN scattering, lies in controlling the required analytic continuation 
of the isoscalar nN amplitude into the unphysical region |86j . which might even be sensitive 
to isospin-breaking corrections [871 [88] . This analytic continuation can be stabilized with the 
help of the low-energy data that have become available in recent years thanks to accurate 
pionic-atom measurements [891 90], leading to a precise extraction of the nN scattering 
lengths HU [92]. A systematic analysis of nN scattering based on this input as well as 
constraints from unitarity, analyticity, and crossing symmetry along the lines of [93H95] will 
help clarify the situation concerning the phenomenological determination of cr^N [96h98] . 

In our numerical analysis, we compare three different methods used to determine the 
scalar couplings and their uncertainties: 

- Method 1: Based on yPT 2 [66], with f^ d determined from (j6j) and from lattice 
QCD. It is well known [99] that the y-nucleon cross section is sensitive to the value 
of / 5 V . In our analysis we adopt the lattice average from [85] : 

/f = 0.043 ±0.011. (13) 


- Method 2: Corresponds to the traditional yPT 3 approach [SB.,[59[[M,.78], where f^ d 
and /‘ ¥ are determined via the three-flavor quantities y and z. In this approach, the 
strange-quark scalar matrix element is defined via 


j,]\f @ttN ITT's 

fs = q - —y ■ 

ZrriN m 


(14) 
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where m s /m = (27.4 ± 0.4) [80], the strangeness content is taken from the relation 
y — 1 — (Jq/ct 7 rTv, with cr 0 = (36 ±7) MeV [ 100] . and z ~ 1.49 is extracted from leading- 
order fits to the baryon mass spectrum ra- This approach introduces uncertainties 
that are difficult to quantify and is sensitive to the precise value of oyjv. The range 
cPriv = (50 ± 15) MeV covering the determinations discussed above [83, 84], [96, [97] 
translates to = 0.2 ± 0.2. Even for moderate values of the cr-term, large values 
R 2 0.25 have been inferred in this way. Such large values are incompatible with 
recent lattice calculations, which provide a more reliable determination of (see (13) 
for a recent average). A determination of the uncertainty bands arising from this 
approach requires us to attach an error to z, which, as argued before, is impossible to 
quantify reliably. Therefore, based on general expectations for the convergence pattern 
in xPT 3 , we simply attach to z a 30% error. 5 

Method 3: Corresponds to the implementation in micrOMEGAs-4. 1.2 |£2] and 
follows the traditional approach in Method 2, 


a function of f. 


but with (14) inverted so that y is 
j v . With lattice QCD input (13) for /; v , this method has reduced 
uncertainties compared with Method 2, but suffers from the fact that f^ d still depend 
on the three-flavor quantities y and z. 


B. Simplified expressions for spin-independent scattering of bino-like dark matter 


Let us now derive analytic expressions for SI y-nucleon scattering in the MSSM. We 
first review the complete expressions due to tree-level Higgs and squark exchange, and then 
simplify them by expanding in powers of u/Msusy- For light third-generation squarks, a 
procedure [102] to extend our results to include the one-loop corrections 121 is discussed 
below. 

The lightest neutra lin o is a linear combination of B , W, and H u ^ interaction eigenstates, 


X = x? = Z*B + Z&W + Z*H d + ZIH U , 
while the neutralino mass matrix is given by 

/Mi 0 -\giv d \giv u \ 
M x = 0 M 2 \g 2 v d -\g 2 v u 

-\giVd \g-2v d o -/i 
V \g\v u -\g 2 v u -/i o / 


(15) 


(16) 


Here M\ ( M 2 ) are the soft SUSY-breaking masses of the bino (wino), /i is the Higgsino mass 
parameter, and v U) d are the two Higgs H u d vacuum expectation values, whose ratio v u /v d is 
denoted by tan (3. Note that while v U)d can be rendered real and positive by an appropriate 
phase redefinition of the Higgs fields, M\ and M 2 are in general complex if the gluino mass is 


5 This is consistent with an analysis m of the quark-mass dependence of octet baryons. 
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X 



FIG. 2: Tree-level MSSM graphs which contribute to the SI cross section for \-quark scattering. 


assumed to be real (as is standard convention). The neutralino mixing in (15) is determined 
by the unitary matrix Z x which diagonalizes M x 


Z}, T Mp P Zj, j — 5jjm^o . 


(17) 


In the squark sector, the squared masses m| s are eigenvalues of the 6x6 matrices in 


flavor/chirality space, 


2 _ f V ] m 2 Q V + vlY u Yl + g UL m 2 z c 2 p —v u (Y u A u + Y u g cot 0) 


Mi = 


-^(Atyt+Y+^cotfl 


rn 


u 


T + gu R m 2 z c 2 0 


Mj = 

a 


_ ( mg + v 2 jY d Y\ + g dL m 2 z c 2f} -v d (Y d A d + Y d fit&n/3) 


-Vd( AjY,J + Y>* tan 0) m 2 D + d + g dR m 2 z c 2 p 


(18) 


Here, the soft SUSY-breaking squark masses are mg, nif/, and m 0 , Y Utd are complex Yukawa 
matrices, V is the CKM matrix, and we have assumed flavor universality & q = Y q A q for the 
trilinear H-terms. We also use s and c for sine and cosine, so that sp = sin/3, c 2 p = cos2/3, 
etc. The weak neutral-current couplings 


9q 


= /.? - e n s 


2 

q?W 


(19) 


are defined in terms of the third component of weak isospin J|, electric charge e q , and 
s 2 v ~ 0.2231. A unitary transformation 


Z%(M 


q)s't'^t't — m 


qjst ) 


( 20 ) 


gives the physical basis with diagonal squark mass matrices, where we adopt the convention 


to order the states in increasing mass. We have also defined the super-CKM basis in (18) 
as the one with diagonal (and in general, complex) Yukawa couplings Y q . [ 104j . 

We have now introduced the necessary ingredients to discuss y-quark scattering in the 
MSSM. The tree-level contributions to the Wilson coefficients C qi are the diagrams shown 
in Fig. [2] For CP-conserving neutralino interactions, these amplitudes were calculated long 
ago in pil - 135] . and extended in [105] to include UP-violating effects. 
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In our conventions, the contributions to C qi due to squark exchange in the s- and u- 
channels at zero momentum transfer read 6 


= l £ 

S= 1 


+ 


_ (m x + m qi ) :2 - m\ + ie (m x - m qi ) 2 -m} +ie_ 


Re {r^*r^ s }, (21) 


where there is no sum over i, and a pole mass m qi enters in the squark propagator. Since 
we work at Mgusy, the running quark masses fh qi must also be evaluated at this scale. For 
the Higgs-exchange contribution we have 7 


*«<%“ = lE Jr Re f r 


k=1 


H, 


Hk \Re lV Hk \ 

XX J t qtqii t 


( 22 ) 


where H\ = H and H 2 = h. We assume that rrih ~ 125 GeV is the mass of the Higgs-like 
resonance found at the LHC pLQ6. FlTT7j and rrih < m h■ In general, the XQiQs and H^XX cou- 
plings appearing in (21) and (22) are complicated expressions involving the mixing matrices 
Z x and ZL Thus the SI cross section is typically determined numerically. However, it is 
known [10811110] that one can obtain analytic results by diagonalizing M x perturbatively 8 
in powers of v/M sUSY . For complex M 12 and / 1 , one hnds to leading order 


zf, = + 0(t. 2 /M| USY ), 

Z21 — 0 (v 2 /M| usy ) , 

— 2 

Z31 — ~ ^2 \m 1 | 2 ^ \fi\ 2 ^ lC P + ^* S ^ + ^(^V-^susy) j 

Z41 = ^2 |Mi | 2 — |/r | 2 ^ lS/3 ^(u 2 /TL|usy) j ( 23 ) 


where 0Mi is the phase of Mi. Note that the presence of a pole at \M\\ = \n\ has no 
physical meaning: it is a consequence of the fact that we assume a bino-like LSP and used 
non-degenerate perturbation theory to diagonalize the neutra lin o mass matrix. 

In Appendix [B] we show how (23) can be used to simplify (21) and (22) if flavor-violating 
effects are neglected, 9 while allowing for non-universal A-terms and squark masses. The 


6 We have i = 1, 2,3 for generation indices and s = 1,..., 6 for squark mass eigenstates. To recover the 
expressions in m, one needs to make the identification Cf JA v>- a 3l etc. 

7 In principle, the CP-odd Higgs can contribute to SI y-quark scattering if ^ or the Yukawa couplings Y q . 
are allowed to be complex. We exclude this possibility in our numerical analysis since we take real ^ and 
the real part of Y q . after threshold corrections are included. 

8 For real M-\ -i and /i, the exact diagonalization of M x is known llllKll3j . although the resulting formulae 
are not simple. 

9 Flavor violation in DM direct detection is strongly suppressed since the effect can only enter via double 
flavor changes / —> j —> f which are experimentally known to be small. Furthermore, the effect of flavor 
off-diagonal entries can be largely absorbed by a change of the physical squark masses. 
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resulting expressions read 


a = i r 


lx r+ >?+ I (T+ If?+1 

g^Ui-^Ui-^Ui + £ A/f 2 _ ,,2 ^ Ui) 


rn _ _d 

°di - o 


54 


C'"’" = 


6 M 2 - \x 2 

-Y r + d+ , I Mi +//tanff 
9 dl di + 6 M ?-^ 2 


(-Mi + /i cot /3) ( —A- + 


+ (L+,i?+) o (L~,R~), 


( L t + *K) 


+ (L + ,R + ) <-> (L-,R~), 


4 Ml - /. i 2 _ 
1 


m l m H 


- (Ml cot f3 + yU)s Q C Q ( - 2 _ 


m 2 m\ 


(Mi + /itan/3) I —% + 


^h,H _ 54 

di 4 Ml - /i 2 _ 


where the squark mixing is defined as 

. = A" + /i cot /3 

while the squark propagators are 

si = 


m\ m 2 H 


— (Mi tan (3 + n)s a c a ( —^ — 






(m x ± m qi ) 2 — m 2 g + ie 


and X fi . = A)] + // tan (3 , 


for 5 = L or i?, 


(24) 

(25) 

? 

(26) 

l] , 

(27) 

(28) 

(29) 


and and are the upper and lower diagonal components of the squark (mass) 2 ma¬ 


trices in (18). In deriving (26 27), we have imposed CP conservation so that the neutralino 


mass parameters M 12 and ft are real. We also take 0 < (3 < n/2 and M 12 > 0 so that 
both signs of /i are allowed. Expressions for C gi in the C P-violating case are provided in 
Appendix [B| 

Note that: 


1. The simplified expressions in (24 25) are valid provided the squarks are sufficiently 


heavy, i.e. if m 1 + m 1 <C m~. This requirement is not met for light third-generation 


squarks, and thus (24 25) must be corrected to account for the one-loop result 


To do so, we follow the prescription adopted in [ IU2 J and replace all tree-level squark 
propagators 

St -> ~ K (±, m q ., m ?, m x ) 


m q ,m q ), 10 

(30) 


(31) 


K(a, m q , rriq, m x ) = § m q [m q (I i - §m 2 / 3 ) - am x (/ 2 - \h ~ § m x h)\ , 

whose form is given in Appendix B of m In the heavy squark limit, the function K 
agrees with to leading order in m- 2 . 


10 


The term proportional to I 3 in (A5) of 1 102) is missing a factor of m q . 
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2. We have made use of the tree-level relation Y q . = m q ./v q in order to obtain (24 27). 
For down quarks, however, this relation can be modified by one-loop graphs which 
induce an effective coupling between d* and the neutral component of H u . These 
corrections [ ITMlllOl 11141 - 1116] are non-decoupling and enhanced by a factor of tan (3 . 11 
For example, the gluino contribution at one-loop modifies the tree-level relation so that 


W = 


m di 


v d ( 1 + e* tan/3) 
where e* ~ -^-mgH*C 0 (mj, m 2 - L , m| fl ) and 

b 2 


(32) 


C 0 (a 2 ,b 2 ,c 2 ) = 


b 2 )(c 2 


b 2) l0g b 2 + 


c 2 )(b 2 


l a 

log — • 


(33) 


Since m qi C qi is proportional to Y qi v q , we can account for (32) by a simple rescaling of 


the Wilson coefficients 


Cl 


,H 


-A 


VdYd, 


■i f^q,H 
U di 


(34) 


where we include corrections [ 109 . 110l II 1414116] beyond the gluino loop (32). These 
threshold corrections feature in our analysis of heavy Higgs H and sbottom contribu¬ 
tions (Sec. Ill) to the SI amplitude. Note that corrections to the light Higgs coupling 
hdd cancel in the relation rh Ui = Y Ui v u . 


III. SIMPLIFIED MODELS: BLIND SPOTS AND ISOSPIN VIOLATION 


We now apply our analytic results (24 27) to four simplified models; each motivated by 
the following experimental and naturalness considerations. Firstly, the ATLAS [TUB] and 
CMS [ 107 j experiments at the LHC have discovered a Higgs boson with SM-like properties 
and a mass below the upper bound < 135 GeV of the MSSM. Secondly, a natural resolution 
of the gauge hierarchy problem requires several conditions [46l [52, [53j to be met: 


• In order to cancel the top-quark correction to the Higgs mass parameter m# , top 
squarks must be light with masses in the sub-TeV range; 

• The gluino mass must be around a TeV in order to prevent radiative corrections driving 
the stop masses too heavy; 

• Light Higgsinos must be present in the spectrum so that tree-level electroweak sym¬ 
metry breaking implies that fi ~ v is satisfied. 


11 In principle, large A-terms can also change the values of Y qi significantly [1101 1117111201 . However, this 
effect drops out in the Higgs-quark-quark couplings where the effective (physical) mass enters. Further¬ 
more, since we assume flavor-universal A-terms in our numerical analysis, the effect cannot be very large 
without violating vacuum stability bounds. 
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FIG. 3: Spectra of the simplified models (A-D) considered in this work. For each model, the LSP 
X is assumed to be bino-like and may be accompanied by a nearly degenerate tau slepton in order 
to produce the observed DM relic density. The SM-like Higgs is denoted by h, while all other 
states are assumed to lie below 1 TeV, including Higgsinos (not shown). From left-to-right, the 
spectra become increasingly more natural as one includes the additional Higgs states H, A and 
third-generation squarks t±,t 2 ,bL- In general, x may be heavier than h, H, and A. 


It has also been observed H2H that naturalness constrains the additional Higgs bosons 
H,H ± ,A to not be too heavy. Barring the gluino, the current experimental bounds on 
the masses of the above particles are rather weak. In contrast, the mass of the gluino 
and squarks of the first two generations are constrained to lie above 1 TeV. Therefore, 
naturalness prompts us to consider the simplified models shown in Fig. [3j where we start 
from a minimal, light particle spectrum necessary to have bino-like DM scattering [model 
(A)] and successively include as active degrees of freedom those particles which are (a) 
required to be light by naturalness, and (b) relevant for DM direct detection. Note that 
due to SU(2)l invariance, the models (C-D) involving two light stops always require a light 
sbottom in the spectrum. (Only if there is a single, mostly right-handed stop, can sbottoms 
be decoupled.) 

In general, a bino-like LSP produces a DM relic density that is too large in most of 
the parameter space considered in Sec. Ill However, the overproduction of bino-like DM 


in the MSSM can be diluted by either s-channel resonance exchange involving Z, h, H , A, 
or x~f co-annihilation with a sfermion / that is nearly degenerate in mass with \. 12 Both 
mechanisms m increase the annihilation cross section before thermal freeze-out and can 
produce the observed relic abundance. In each of the models shown in Fig. [3j the relic 
density constraint may be satisfied by either mechanism or, if necessary, by extending the 
spectrum to include a tau slepton t which generates additional co-annihilations mm- 
Since the f mass can be tuned without affecting naturalness or DM direct detection, we do 
not consider the DM relic density constraint in our subsequent analysis. 

Similarly, we do not consider the constraint from the anomalous magnetic moment of 
the muon a^, whose world average is dominated by the Brookhaven measurement m- 
The resulting value deviates from the SM prediction by 3-4cr, depending on the details of 


See e.g. DU for a detailed analysis of these effects in the pMSSM. 
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the evaluation of the hadronic contributions [12611127] . Recent developments in the eval¬ 
uation of the SM prediction include: the QED calculation has been carried out at 5-loop 
accuracy cm after the Higgs discovery [1061 107 ] the electroweak contribution is com¬ 
plete at two-loop order [129] 130] . and hadronic corrections have been considered at third 
order in the fine-structure constant [ 1311 . 132] . Although an improved determination of 
the leading hadronic contribution, hadronic vacuum polarization, mainly requires improved 
data input, see [ 1261 I133H135] , the uncertainties in the subleading hadronic-light-by-light 
contribution have been notoriously difficult to estimate due to substantial model depen¬ 
dence [i26i Tm nasi- Recently, data-driven techniques have been put forward to reduce 
the model dependence based on dispersion relations and a first lattice calculation 

has become available [141]. All these efforts are motivated by two new experiments, at 
FNAL [T42] and J-PARC [143] . which each aim at improving the measurement by a factor 
of 4 and thus help clarify the origin of the discrepancy between experiment and the SM 
prediction. 

Should the discrepancy persist, an explanation within the MSSM is possible provided 
certain assumptions are made about the SUSY parameters entering the smuon, chargino, 
and neutralino mass matrices. If these parameters are all equal to Msusy, then a positive 
contribution to a p requires sign (/iM 2 ) > 0 since the dominant one-loop amplitude scales 
approximately with /iA7 2 tan [3 /M| usy ; see e.g. the review [144] and references therein. In 
the blind spot regions where /i < 0, this condition would require us to relax the assumption 
that M 2 > 0. However, the requirement sign (/iM 2 ) > 0 does not necessarily apply if the 
SUSY mass parameters are non-degenerate. For example, it has been shown |145fH47] that 
a positive contribution to a p can arise if \Mi\,rriji R <C |M 2 |,m^ i , in which case /i and M 1)2 
must have opposite sign. The key point is that neither the sign of M 2 nor the smuon 
masses are relevant for our analysis of SI scattering, so it would be possible to account for 
the experimental value of a p by a suitable choice of these parameters. Furthermore, the 
discrepancy could also be explained by large A p terms [1181 114811149] not correlated with 
DM scattering. 

We conclude this section by anticipating a key result of our analysis: isospin-violating 
effects can be magnified in the proximity of blind spots, where the SI direct-detection cross 
section lies below the lower bounds set by the irreducible neutrino background. For these 
parameter-space configurations, the SI amplitude itself becomes tiny and hence more suscep¬ 
tible to small variations in the input quantities, such as corrections from isospin breaking. 
In particular, the ratio of proton and neutron SI cross sections becomes very sensitive to the 
values of the scalar matrix elements and their uncertainties 5f n%p , 



so that the overall uncertainty on f n /f p can become large near blind spots where f p ~ 0. 
In each of the four simplified models (A-D), we examine the amount of isospin violation 
associated with the three methods of Sec. Ill Al 
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tan/3=2 



tan/?=10 



tan/3=5 



tan/3=50 



FIG. 4: Current and projected limits on SI y—xenon scattering due to h exchange with tan (3 = 2 
(top left), 5 (top right), 10 (bottom left), and 50 (bottom right). The pink band shows the existing 
constraints from LUX ED, while projected limits from XENON1T |22| and LZ [23] are given by 
the blue and orange regions respectively. The blind spot where the SI cross section vanishes is 
denoted by the red line and lies within the irreducible neutrino background (r'Bc) shown in gray. 
We assume the LSP is bino-like (Fig. [3]), so do not consider the triangular, hatched region where 
fjL < M\ and \ becomes Higgsino-like. 


A. SM-like Higgs exchange 

We begin by considering the minimal particle spectrum for which an observable SI cross 
section is possible. From the rightmost diagram in Fig. [2j it seems reasonable to conclude 
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that the SM-like Higgs and the bino-like LSP is sufficient in this case. However, in the limit 


rriA rriz, (26) and (27) become 


Ct = C h d . = 

U>i LLi 


9i 


4m); Mf — fi- 


-(Mi + fj,s 2 p) , u,i = u,c,t, di = d,s,b, 


(36) 


and thus the scattering amplitude decouples with the Higgsino mass //. It follows that a 
measurable cross section due to Higgs exchange implies the presence of light Higgsinos in 
the spectrum, thereby satisfying one of the minimal naturalness requirements. Although 
this feature does not prevent the reintroduction of fine-tuning in the MSSM altogether, it 
becomes relevant in our subsequent analysis where light stops are added to the spectrum. 


To compare (36) to data, we first note that Cg. vanishes when 


Mi + /US2/3 — 0 , 


(37) 


and thus a blind spot arises in the SI cross section provided // is negative. The prospects for 
constraining this feature (37) have been extensively analyzed [SU] for y-nucleon scattering. 


To examine isospin violation, however, we need limits on x~ nu cleus cross sections, so we use 
0 in order to constrain the relevant parameter space. 

Let us first consider the limits associated with ( |36| when the scalar matrix elements fg 
of Method 1 are employed. In Fig. [4j we update the results from [30J and show constraints 
for various values of tan/3 in the (/qMi) plane from current and upcoming xenon experi¬ 
ments. For yU > 0, we fold that only a narrow strip is excluded by the existing limits from 
LUX [21J, while the projected reach from XENON1T [22] and LZ [23] will probe most of the 
naturalness-preferred region where // is of order v. As tan f3 is increased, the term oc /i in Cg. 
is suppressed, thereby weakening the direct detection limits. If no signal is seen at LZ, then 
the allowed parameter space is focused towards tan /3 = 50 and values of Mi < 200 GeV. In 
the n < 0 region and for small tan/3, the naturalness-preferred values of fi occur at |/i| ~ Mi 
and are concentrated near the blind spot. Although the irreducible neutrino background 
make this region difficult to probe experimentally, larger values of tan /3 decrease the blind 
spot slope, so that natural values of /i become allowed for Mi < 300-400 GeV. 

By taking a slice through the (/i, Mi) plane, we can also extract the limits due to a small 
mass splitting /i| — Mi = 80 GeV between the bino and Higgsinos. This choice is motivated 
by the current CMS results (150) on same-flavor opposite-sign dilepton searches. Here CMS 
sees a 2.6cr deviation which can be explained by a heavier neutralino decaying to a lighter 
one. Fig. [5] shows the resulting constraints, where we plot the SI cross sections as a function 
of the bino mass. For /i > 0, the limits from LUX are stringent, with values below Mi ~ 600 


GeV excluded. The strength of these limits is due to an enhancement in the amplitude (36) 


from both a nearly degenerate denominator and lack of interference in the numerator terms. 
For jjt < 0, there are no constraints from LUX, although XENON1T and LZ will exclude 
the whole parameter space in the absence of a DM signal. 

We now examine the hadronic uncertainties associated with each of the three methods 
discussed in Sec. IIA For h exchange, the Wilson coefficient (36) is independent of quark 


16 












tan/?=10, \iJ\-M- i =80 GeV 



FIG. 5: SI y-xenon cross sections for h exchange with tan/3 = 10 and a small mass difference 
|/x| — M\ = 80 GeV between the Higgsino and bino. The solid (dashed) black line corresponds to 
// > 0 (/r < 0). Shown are the limits from LUX [21], XENON1T [22], and LZ [23] , with the same 
color coding as in Fig. [4] 



FIG. 6: Amount of isospin violation in terms of f n /fp due to h exchange in \-nucleus scattering. 
The colored bands correspond to the la uncertainties associated with the different determinations 
of the scalar matrix elements discussed in Sec. IIA The blue band corresponds to Method 1, 


while the orange and green bands correspond to Method 2 and Method 3 respectively. 


flavor, so the SI amplitude ([2]) factorizes 



Evidently, the resulting SI cross section is sensitive to the value of /'/ v , with a dramatic 
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effect observed [99] on the regions of excluded parameter space when the typically large 
value ~ 0.25 of Method 3 is replaced with much smaller determinations (13) from the 


lattice. We emphasize that this sensitivity is also present in any analysis of isospin violation, 


where f n / f P is the quantity of interest. For the present discussion, (38) implies that the ratio 


fn 

f P 


m r 


2 + 7 Zfo 


m p j 2 + 7 Y,fi 


(39) 


is independent of C ^., and thus isospin violation is entirely determined by hadronic quantities. 


In Fig. [6] we compare the uncertainties on f n jf v as a function of a nN . For Methods 1 and 
3, we find stability across a large range of cr nN values, with isospin violation allowed at 
around the five and ten percent level respectively. As noted in [T5T] , this stability is due 


to the fact that the constant term of | in (38) dominates the remainder whenever ff r is 
fixed by lattice input. In contrast, the yPT 3 formalism of Method 2 produces a strong 
dependence of f!f on c+jv, which in turn affects f n jf P - From Fig. [6j isospin violation greater 
than 50% is allowed, in marked contrast to the precision of Method 1. This example clearly 
demonstrates the huge uncertainties associated with Method 2, which, however, is still used 
in the literature [ TOl I152j . 


B. Light and heavy Higgs exchange 


m 


Let us now extend model (A) to include the heavy Higgs bosons H, A, H ± [model (B) 
Fig. 1- The inclusion of these additional degrees of freedom is motivated by natural¬ 
ness [121], however, only H contributes to the SI cross section (Fig. [2]). 

From our simplified expressions (26 27), we see that the couplings to up and down quarks 


differ by a factor of tan/3, but are identical 13 among different generations i = 1,2,3. As a 
result, the SI amplitude may be expressed as 


In 

m N 


= ct H u, 


N 


+ C h y H D 


N 


(40) 


where 


Ov = /" + 2/£' and D N = ff + ff + f Q (41) 

collect the scalar coefficients associated with the up- and down-type Wilson coefficients. A 
blind spot occurs if the condition 

C h ^U N + C h d fD N ~ 0 (42) 

is satisfied, and the resulting suppression of the SI cross section has been identified nu¬ 
merically j39[ [55Tf5T] and further studied analytically [38] . In the latter case, an explicit 


13 


Up to threshold corrections (34), which enhance C^ by tens of percent at large tan/3, 
does not have a large impact on the numerical analysis. 


Their inclusion 
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FIG. 7: Dependence on a n N i n the ratio of the hadronic terms Un and Djy pre-multiplying the 
up- and down-type Wilson coefficient for h,H exchange (40). Shown is the case for N = neutron, 
with similar results obtained for N = proton. Color coding as in Fig. [6} 


formula [38] for the blind spot can be found for moderate to large values of tan f3 and 
rriA > m h : 

9 1 

~ 0. (43) 


2 1 
— (Mi + /xs 2| a) + ptan 


m 


m 


H 


In effect, (42) has been recast as an interference condition between the h and H amplitudes; 
a feature which has important consequences for isospin violation in the MSSM. As with h 
exchange, negative values of fi are required in order to generate the blind spot. However, 
note that in the vicinity of (37), the first term in (43) is suppressed, so in some cases the 
contribution from H exchange may dominate the scattering amplitude [38 ]. 

A crucial step in the derivation [38] of (43) is the observation that Un ~ Dn numerically. 
Deviations of D^/U^ from unity have the effect of shifting the location of the blind spot 
(43), so it is necessary to determine this ratio precisely. In Fig. [Tj we display the sensitivity 
of D n /Un to cpriv for each of the three methods of Sec. II A| Similar to our analysis of h 
exchange (Fig. [6j) , we find that Methods 1 and 3 are stable across a large range of a n N values, 
with Dn/Un ~ 1 tightly constrained. In contrast, Method 2 exhibits a strong dependence 
on a wN and for a nN > 45 MeV, the location of the blind spot (43) can get shifted by a factor 
of eight or more. These findings illustrate again the importance of using a well-controlled 
framework for the hadronic input quantities. 

Let us now examine the experimental limits associated with y-xenon scattering. In 
Fig. [8j we show constraints in the (tua, tan/3) plane for two benchmark values of Mi and 
/i. We find that as the mass splitting between Mi and /i is decreased, the limits become 
significantly stronger. This is because the amplitude C^ H scales like ~ l/(Mf — p 2 ), so the 
naturalness requirement of light Higgsinos implies strong constraints on the SI cross section. 
We also find blind spots similar to those previously identified [38], [SSj, I55H57] . and see that 
the strongest limits are due to H, A —> t + t~ searches [E3] as one approaches ([43]) from 
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Mi =220 GeV, p=-440 GeV /W 1 =220 GeV, M-M, =80 GeV 




m A [GeV] 


m A [GeV] 


FIG. 8 : Current and projected limits on SI y—xenon scattering due to h, H exchange with different 
benchmark values for Mi and p. Excluded regions and the blind spot are color-coded as in Fig. [4j 
with the cross-hatched region in dark-blue corresponding to CMS limits [ 153] on -> t + t~. 
The region to the left of the dark-red dashed line at mA — m#+ — 480 GeV is excluded by 

B^x sl PB]. 


below. 

What about isospin violation in this model? Unlike single h exchange, where f n /fp is 
entirely fixed (39) by hadronic quantities, the blind spot (43) for light and heavy Higgs 
bosons involves destructive interference between the respective amplitudes. In general, we 
find that isospin violation can be enhanced in the vicinity of such blind spots because f n /f p 
becomes sensitive to the scalar matrix elements and their uncertainties (35). This is evident 
in Fig. [9j where the central value of f n /f P (determined by Method 1) reaches ~ 15% as the 
blind spot is approached with increasing mA- 

In Fig. [lO] we compare the amount of isospin violation allowed by Methods 1-3. As 


observed in single h exchange, the uncertainties associated with Method 2 are large, and 
differ by a factor of two or more for mA ~ 400 GeV. For rriA ^ 400 GeV, a comparison 
between the Methods is obscured by the fact that the location of the blind spot is shifted 
depending on deviations from Dn — Un (Fig. [7]). 

In Fig. 11 we display the allowed ranges of isospin violation due to Method 1 for two 
values of tan (3. For tan f3 = 10, the recent limit from B —>■ X s ~f [ 154] excludes the region 
below the blind spot at ~ 500 GeV. Above the blind spot, the absence of a signal at 
LZ would imply that isospin violation as large as 10% becomes allowed. For tan (3 = 20, 
the current limit from LUX [21] allows around 10% isospin violation, although this occurs 
at a value of rriA ~ 200 GeV already excluded by the limits from H, A —>• r + r~ [153] and 
B —>■ AT ,7 [154] , The absence of a signal at XENON1T would allow ~ 20%, while at LZ 
this would imply that isospin violation as large as 40% is allowed within the uncertainties as 


20 





























m A [GeV] 


FIG. 9: Pseudoscalar Higgs mass m A dependence of the SI y-xenon cross section (black) and the 
central value of f n /fp (red) as determined by Method 1. 


M-j =220 MeV, ^=-440 MeV, tan/3=10 



m A [GeV] 


FIG. 10: Amount of isospin violation in terms of f n /fp due to h, H exchange in y-xenon scattering. 
The shaded regions show the uncertainty on f n / f p due to each determination of the scalar matrix 
elements listed in Sec. IIA Color coding as in Fig. [6j 


one approaches the blind spot from below. As illustrated in Fig. [9j it is important to note 
that not only the allowed range, but also the central value of f n jf p can increase as the blind 
spot is approached. In consequence, the absence of signals in SI DM searches pushes the 
parameter space into blind spots, at which f n jf v may become large and thus the accurate 
determination of 5f n and 5f p becomes paramount. A comparison for other nuclear targets 
can be inferred by taking the limits on f n jf p and comparing against Fig. [lj 
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M, =220 MeV, ju=-440 MeV, tan/3=20 


Mi =220 MeV, p=-440 MeV, tan/?=10 




m A [GeV] m/i [GeV] 

FIG. 11: Current and projected results for the measure f n /f P of isospin violation arising from h, H 
exchange in SI y-xenon scattering for tan (3 = 10 (left) and 20 (right). The dark blue region denotes 
the uncertainty on f n /f p as determined by Method 1. The existing and projected experimental 
limits are color coded as in Fig. [4j with the irreducible neutrino background (z^BG) shown by the 
central gray band. The region to the left of the dark red dashed line at mA — 480 GeV is excluded 
by B Xs'y |154j . 


C. SM-like Higgs and light squark exchange 


In the previous subsection, we investigated parameter configurations (42) where the h, H 


amplitudes interfere destructively and observed that isospin violation can be enhanced in 
the proximity of these blind spots (Fig. 0- Next, we examine if blind spots still exist once 
third-generation squarks are added to the spectrum of model (A). 

The effects on the SI amplitude due to squarks from the first two generations were consid¬ 
ered in [30] (including h exchange) and shown to be small due to the stringent limits from 
LHC searches. However, the existing limits on third-generation squark masses are much 
weaker, so that effects from stops and sbottoms can be significantly larger. 

The simplest model, i.e. with minimal particle content, would involve a single, mostly 
right-handed stop t r. However, as one can see from (24), this contribution is not tan (3 


enhanced and thus the h contribution (36) dominates the SI cross section. Therefore, we 


consider a spectrum where h and £ 1j2 are the dynamical degrees of freedom [model (C) in 
Fig. [1]. Since a left-handed stop is always associated with a left-handed sbottom b L , the 
sbottom contribution must be included as well. Although this does not increase the number 


of free parameters, we can see from (25) that the sbottom amplitude is tan [3 enhanced and, 


crucially, can compete with the Higgs contribution to the SI amplitude. Note that while 
the Higgs amplitude vanishes with decoupling Higgsinos (thereby violating the minimal 
naturalness conditions), this is not the case for the sbottom contribution, which possesses a 
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term proportional to fi. 

Using our simplified expressions, we find that a blind spot occurs if the condition 


6 

m 2 h 


(Mi + /1S2/3) 


EU + 3/q 


Ml 


m 


-(Mi + fitanf3)f Q ~0 


bL 


(44) 


is satisfied. Here, we have ignored the numerically small stop contribution C f t and ap¬ 
proximated the effects due to sbottom loops (31) by the tree-level propagator. This latter 


approximation is for illustrative purposes only, and in our numerical analysis we use the 
exact one-loop expressions. Using the scalar matrix elements from Method 1, we have 
Y q fq + 3 fg ~ | and fq ~ A. Therefore, for moderate to large values of tan/3 the blind 
spot condition simplifies to 


30 , „, 

—+ /is 20 ) + /itan/3^ 
mi m 


b L 


-M 2 


~ 0. 


(45) 


As expected, this blind spot shares common features with the one found |38j for h , H ex¬ 


change (43): it requires negative values of /i, so that the couplings to h are suppressed and 


destructive interference between the h and bL amplitudes can occur. However, larger values 
of |/x| are required in order to overcome the factor of 30 in the Higgs amplitude. 

Before determining the experimental limits on this model, let us consider the size of the 
parameter space. The h amplitude depends on M 1; fi, and tan/3, so we need to add the 
parameters (mg) 33 , (m^) 33 , and X t of the stop mass matrix. As noted above, the sbottom 
contribution does not involve additional parameters since the left-handed sbottom mass is 
given by 

m l L = ( m Q)s3 + m b ~ (| - \swWzC 2f} . (46) 

To reduce the number of free parameters we fix (m<g) 33 « (m{/) 33 , in which case the left- and 
right-handed entries in the stop mass matrix become nearly degenerate, while the physical 
mass eigenvalues read 

m~ i >fa = m~ L =F m t X t . (47) 


This allows us to express our simplified expressions (24 25) in terms of the physical masses 


and compare with collider limits in the (my, Mi) plane. In order for light stops to generate 
the correct Higgs mass in the MSSM, we assume these states are mixed in such a way so 
as to give a maximal contribution to the Higgs mass. This is achieved by noting that the 
one-loop stop contribution to the Higgs mass 


m h 


2 2 1 

m Z c 2/3 + 


3 m 
47T 2 V 2 


4 r 


, m~ X 2 
hr H—ry 

mr 


1 - 


12 m 2 - 


(48) 


is maximized at \X t \ = X) nax = \/6 my, where m 2 = mymy 2 is the average stop mass. 
Since m/> is bounded at tree-level by m/, requiring m^ ~ 125 GeV implies a lower bound 
m 2 > 550 GeV for X 4 max . While this is the main source of fine-tuning in the MSSM, it can 
be easily evaded in e.g. non-minimal SUSY models, where the correct Higgs mass can be 
obtained via non-decoupling D-ternrs [155] . or in the next-to-minimal supersymmetric SM 
with special parameter choices [156] . 
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FIG. 12: SI y-xenon cross sections as a function of the lightest stop mass rn^. The solid lines 
correspond to different values of the Higgsino mass //, while the pink dashed line is the existing 
limit from LUX. The blind spot at rri^ ~ 160 GeV is shown in red. 


AT,=100 GeV, tan/3=10 



In Fig. 12 the SI cross section is displayed as a function of the lightest stop mass mi 
for tan /3 = 10 and several values of /i. We find that a positive value of /i is excluded by 
LUX, while negative values become progressively harder to constrain as the mass difference 
between the bino and Higgsinos is increased. The blind spot is clearly seen for fi = —4 Mi 
and occurs at a light stop mass rri^ ~ 160 GeV. 

In Fig. [13] we show the interplay between collider and DM direct-detection limits in 
the plane for two values of tan/3 and /i. We find that in the absence of blind 

spots (/i = —2 Mi), LUX excludes the M { < 50 GeV region across a large range of stop 
masses. These limits will be significantly improved if XENON1T and LZ do not detect 
a DM signal, with whole regions below M\ « 300 GeV and 500 GeV excluded by the 
respective experiments. In these cases, the direct-detection limits surpass those derived 
from the ATLAS searches for stops and sbottoms. 

In the blind spot region (/i = —4Mi), the DM limits are considerably weakened, with 
LZ excluding most of the M\ 50 MeV and 150 MeV regions for tan/3 = 10 and 25 
respectively. As noted in Sec. Ill B isospin violation can be enhanced in the proximity of 
blind spots which arise from destructive interference in the amplitude. Since this is a generic 
feature of SI scattering, it follows that isospin violation can also be large near the blind spots 
shown in Fig. 13 Although these regions are excluded by collider limits, we have not ruled 
out the possibility that blind spots for this model occur in viable regions of parameter space. 
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FIG. 13: Current and projected limits due to light Higgs h and squark exchange in x~ 

xenon scattering with tan (3 = 10 (top row) and tan (3 = 25 (bottom row) and two benchmarks for 
negative p. The figures show the DM constraints from LUX [21] and XENON1T [22] (color coded 
as in Fig. [4]), while limits from direct searches for stops and sbottoms at ATLAS |157j are shown 
in green. The blind spot is shown in red and lies within the neutrino background (i'bg) shown in 
gray. The hatched region corresponds to the case where t\ becomes the LSP. 


D. Generic Higgs and light squark exchange 

In this Section, we consider the effect of adding H to the particle content, so that the 
active degrees of freedom are h,H and £ 1i2 ,&l [model (D) in Fig. [ 3 ). This is the most 
“natural” model studied in this article since a large value of m H would also require fine 
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tuning [121]. 

To derive an analytic formula for the blind spot, we follow the same steps used to obtain 
(44). For moderate to large values of rriA > rrih and tan/3, we find that a blind spot occurs 

15 


whenever 


30 


—2 {Mi + P-S2/3) + /U tan /3 I — ^ —t~ 


m 


m 


H 


m~ 


b L 


-M\ 


~ 0 


(49) 


is satisfied. In this case, the inclusion of H has the effect of shifting the location of the 
blind spot found for h and squark exchange (45). In particular, negative values of // are still 
required. 

To examine the limits on this model, we must also consider flavor observables since H and 
light stops contribute to B —* X s j and B s —> n + n~. To evaluate the flavor constraints we 
use SUSY_ flavor- 2.51 [15811160] and implement the NNLO SM calculation by constraining 
the ratio 


-Rsusy {B -A- X s q) = 
to lie within the allowed range for 


Br (B —>■ X s7 ) 


MSSM 

SUSY FLAVOR 


\SM 


Rexp(B -a- X s7 ) — 

Here we use the recent calculation of [ 154 ] 


Br(B —y X s7 )g USY flavor 
Br(B —>■ A" s7 )exp 


Br (B —>■ X s7 ) SM 


NNLO 


Bt(B X s7 ) SM = (3.36 ± 0.23) x 10 


-4 


(50) 


(51) 


(52) 


to incorporate the NNLO SM prediction. For the experimental value we use the PDG 
average [161] within 2cr uncertainties, 


Bt(B -a- X s7 )exp = (3.41 ± 0.21 ± 0.07) x 10 


-4 


(53) 


We add the theoretical error linearly with the experimental one, so that -Rsusy is required 
to lie within the interval 


Kw{B X. s7 ) < R S \jsy(B -a X s1 ) < R^xp(B -a X s1 ) . 


MAX/ 


(54) 


For B s —> /r + /i we adopt the same procedure to impose limits on the relevant SUSY 
parameter space. Here the SM prediction HS2] is 


Br(B s -A /t + /i')sm = (3.65 ± 0.23) x 10“ 9 
and has to be compared against 

Br (B 8 exp = (3.1 ± 0.7) x 10“ 9 . 


(55) 


(56) 


In Fig. 14 we display limits in the (m 7] , M \) plane for tan/3 = 10 and several values of 
rriA■ Since this choice of tan (3 corresponds to a horizontal slice through the h, H parameter 
space (Fig. [8J, the effect of increasing is to probe the effect of the h,H blind spot (43) 


26 


















A7=-2Mi, tanjS=10, m A =300 GeV v=-2M,, tan/3=10, m A =515 GeV 



m- [GeV] 

‘ 1 


m- [GeV] 

‘1 



100 200 300 400 500 600 700 800 


m- [GeV] 

*1 


FIG. 14: Current and projected limits in the (m^, M \) plane from h,H and ii, 2 )&L exchange in 
X'-xenon scattering. The three plots represent different values of rriA, for fixed tan (5 = 10. For 
4“ax < o, the regions between the purple dashed lines are excluded by B s —> , while regions 

to the left of the dark red dashed lines are excluded by B —> Xg'y. Excluded regions from direct 


detection are color-coded as in Fig. 13 


from below. For tua = 300 GeV and away from the blind spot, we find that LUX excludes 
the band below Mi « 100 GeV. In this case, the limits from B —y X s ~/ and B s —y n + n~ 
provide a complementary and stringent constraint: compatibility with both observables and 
LUX only leaves a small region of parameter space viable. XENON1T and LZ will carve out 
most of the remaining parameter space, providing a very strong constraint on the light 
scenario. However, for heavier values of m^, the constraints from flavor and direct detection 
weaken considerably and here the collider bounds become the dominant constraint. This is 
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FIG. 15: Current and projected limits in the (m^, M \) plane from h,H and ii, 2 )&L exchange in 
y-xenon scattering. In the figures, the value of m,A is increased for fixed tan (3. Excluded regions 
are color-coded as in Fig. [T4j with B —> Xg'y ruling out both values of tua entirely. 


particularly evident in the second plot of Fig. [14j where the blind spot suppresses the SI 
cross section and the moderate value of tiia reduces the tension with B —> X s ~f entirely. 

For iriA = 600 GeV, there is no 


We consider the effect of increasing tan/I in Fig. 15 


blind spot in the physical region, with most of the area below ~ 450 MeV excluded, 
while for = 750 GeV a blind spot does occur. This latter region is allowed by current 
collider limits, but excluded by B —> X s j. However, one should keep in mind that the 
flavor bounds included in Fig. 14 (and Fig. 15) are the least rigorous ones as they can be 
evaded if some of the underlying assumptions are relaxed: in the presence of non-minimal 
sources of flavor violation the bounds can become weaker. In fact, a mass splitting among 
the left-handed squarks deviates from naive minimal flavor violation since there are either 
off-diagonal elements in the up or in the down sector of the squark mass matrix (or in both 
simultaneously). Furthermore, relaxing our assumption that the left-handed bilinear terms 
are diagonal in the down basis (18) would lead to additional effects in flavor observables. 


IV. CONCLUSIONS 


In this paper we examined four simplified models in the framework of the MSSM where 
all but a handful of superpartners are decoupled from the spectrum. We started from the 
minimal model necessary to provide a viable DM candidate and sequentially added particles 
to render the spectrum more natural. The key result, analytic expressions for the Wilson 
coefficients relevant for the Higgs- and squark-exchange contribution to spin-independent 


WIMP-nucleon scattering, is summarized in (24 27) 












As the main application of our scheme, we studied the amount of isospin violation gener¬ 
ated by single-nucleon contributions to the spin-independent WIMP-nucleon cross section. 
In general, isospin violation is a rather small effect: for pure h exchange, the amount of 
isospin violation is « 5%, i.e. in line with common expectations and results in two-Higgs- 
doublet models [163j. Beyond single h exchange, however, the effect can be enhanced in 
the proximity of blind spots. In such cases, the blind spots occur (at a given order in 
perturbation theory) due to destructive interference among different contributions to the 
SI amplitudes. As a consequence, small variations of the amplitude become increasingly 
important. Although the SI cross sections are strongly suppressed in the vicinity of blind 
spots, as direct-detection experiments become increasingly more sensitive, MSSM models 
are pushed towards these corners of parameter space. For instance, we find that for h and 
H exchange, the projected limits from LZ [23] allow isospin violation to be as large as 40%, 
which increases rapidly as one approaches the irreducible neutrino background. In this way, 
precise predictions for isospin violation are essential to relate future direct-detection data to 
MSSM predictions. 

In our simplified models, the source of isospin violation originates purely from the SM; 
it is the blind spots that make these small differences prominent. This situation is unlike 
e.g. in Z' models which introduce new sources of isospin violation beyond the SM. There¬ 
fore, an accurate evaluation of isospin violation requires a careful assessment of the nuclear 
input quantities and the associated hadronic uncertainties. We demonstrated this point by 
comparing different methods to determine the proton and neutron scalar matrix elements 
that are currently used in the literature. While the traditional approach based on yPT 3 
relations suffers from large and, in part, uncontrolled uncertainties, the hadronic input can 
be accurately evaluated by using the two-flavor formalism developed in pB]. In particu¬ 
lar, we showed that in the three-flavor framework, depending on which input is used for 
the strangeness-related quantities, incorrect conclusions concerning both central values and 
uncertainties can occur. 

We also extended our models to include light stops and sbottoms in the spectrum. Again, 
for certain corners of the parameter space, cancellations occur that suppress the amplitude 
and led us to the identification of new blind spots. We identified these blind spots analytically 
in (45) and (49). Furthermore, the interplay between DM, collider, and flavor limits was 
studied, finding that the inclusion of the latter tends to exclude configurations with blind 
spots which are allowed by collider bounds. Only for tan/3 = 10 and near the blind spot 
generated by h and H exchange, did we find a blind spot consistent with all constraints 

(Fig-0- 
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Appendix A: Perturbative diagonalization of the neutralino mass matrix 


In this Appendix, we diagonalize the neutralino mass matrix M x (16) via a pertur¬ 


bative expansion in n/MgusY- Following |108H110] . we first consider the diagonal matrix 
Z X 'M X 'M X Z X , where to leading order in v/Ms\jsy we have 


/ 




1 

0 


0 

1 


g\v (MiC 0 +n*sp) g 2 v (Micp+p*sp) 

V2 72 |M!| 2-| M |2 

g 2 v (MlSp+/J,*Cp) g 2 v (M^Sp+p.* Cp) 

\ 72 |JUi | 2 —|/i | 2 V 2 |Mi| 2 -H 2 


g\v {M^cp+psp) g lV (M^sp+p,cp) \ 

72 \Mi\ 2 ~M a 72 |Mi| 2 -|ai| 2 

_ Q2V (W Cp+P s p) g 2 y (M 2 Sp+pcp) 

72 \Mi \ 2 — \g \ 2 72 IMaP-17 2 

1 0 


0 


(Al) 




Although Z x diagonalizes the square M X ^M X , we need to perform two additional rotations 
in order to make Z xT M X Z X real and diagonal: 


7 X - 
— 



1 

0 

0 

0 

\ 


0 

1 

0 

0 



n 

n 

1 

1 




VJ 

V2 

V2 


V 

0 

0 

1 

V2 

1 

V2 

J 


and 


Z} = 


(J2 0 0 0 \ 

0 e -i<t>M 2 ! 2 o 0 

0 o 6 -#m / 2 o 

V o o o e -^/ 2 ) 


, (A 2 ) 


where <j>Mi,p, is the phase of M 1j2 and 7 respectively. The resulting mixing matrix is given by 


z* = Z*z*zx , 

from which we deduce the relevant components for the lightest neutralino 

zi i = + 0(v 2 /Mg US y) , 

Z 21 — 0 (w 2 /-M| ugY ) , 

e ~2 ‘t’Mi 


(A3) 






9 iv 


y/2 \Mi\ 2 


-5 2 


9iv 


(MiCjg + 9* sp) + 0(v /Af SUSY ), 


:(MiSp + 9 * 0 / 3 ) + 0 (u 2 /M| usy ) • 


a/2 \Mi\ 2 — |//| 2 

Appendix B: Analytic expressions for x ~nucleon scattering 


(A4) 


This Appendix concerns the derivation of the analytic expressions (24 27) for the Wilson 
coefficients C qi appearing in the scalar %-nucleon couplings (J2|. The derivation involves an 
analysis of the spin-independent (SI) amplitude for \Qi XQi scattering due to tree-level 
Higgs and squark exchange (Fig. [2j) . 
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Squark exchange —Consider first the contribution due to squark exchange, where the 
zero-momentum propagator for the s- and u-channels is denoted by 

Dqiqs ± m q .) 2 — m 2 - s + ie ’ 

and we define 

i(Vf s *P R + r f s *P L ) (B2) 

as the Feynman rule for the XQiQs coupling. Then for spinors u, and Vi carrying momentum 
Pi, the s-channel amplitude is 

a« = £ s,[i(rf'J 5 * + rf--p L )]u 2 [t Ad s4t(rf p L + r^p* 

s=l 
. 6 

= -5 E D ti. {n i -'rf-(* i rPLV 3 )(u i%t P R U 2 ) + Tf‘ , Tf-(v 1 YPRVs){u i '1 R P L U 2 ) 
+T q l qa *Y q ^ qs [(v 1 P R v 3 )(u 4 P R u 2 ) + \{vi^v 3 )(u 4 a^P R u 2 )] 

+T q j i qs *T q l qs [(viP l v 3 )(u 4 P l u 2 ) + |(hio rAtI/ u 3 )(M40-^P L M 2 )] } , (B3) 

where the Fierz identities 


( PL, R )ij(PL, R )kl = \{PL, R )il{PL, R )kj + ixvPL, R )kj , 

( PL, R )ij(P R ,L)kl = ^(,'y ll PR,L)il( r )'iiPL, R )kj , (B4) 

have been used to obtain the final equality. Similarly, in the '(/-channel we find 

AS = - E ®3[i(r TPl + rs*-P«)]«2 [*Ad uMpT’Pr + 

S=1 
. 6 

= \ EAs. {rf-rf■ •(v il <‘P R v l )[u ilr p L u 2 ) + r« i -rr(®37'‘Pi't'i)(s47„PR«2) 

Z 5=1 

+Y q i qs Y q ^ q9 * [(u 3 P L V\){u 4 P L u 2 ) + 1(v 3 ^ u Vi)(u 4 (t^P l u 2 )\ 

+r* q T q £ q '* [(v 3 P R Vi)(u 4 P R u 2 ) + 4(v 3 (T M ‘ / 'Ui)('i/4(T Mi ,P R U2)] | , (B5) 

so neglecting the spin-dependent terms involving 7^ and 7^75 gives, 

. 6 

^Uls. -jE [As. + Ad R* {rr-rr**} («3®i)(“4« 2 ) + 75 terms. (B6) 

^ 5=1 


At zero-momentum transfer, the terms involving 75 are suppressed and so we can read off 
the Wilson coefficient for the operator xXQiQf 
6 


™«.A = l E Ah. + Ad Re {rfTS 


<?S* 


S=1 

1 6 

ly 

8 ^ L ( m X + ™?i) 2 - m i + ^ 1 (™X - - m i + ie . 


Re <! rjf J T| 9a * 


(B7) 
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where there is no sum over q t and fh qi is the running quark mass. In the literature, the quark 
mass in D^- s is often neglected, in which case the s- and ^-channel amplitudes coincide with 
each other. Note that by substituting the couplings [103J 


dj. ds 


^did s 

R 


d* 
is 


= V2g 2 {\Z^ - |tan 0 w Zj C 1 )Z s 
= - ^92 tan e w Z*’Z?;, s - Y dl zgzf ;, 


\z* yX yd* 

* di^ 31^2+3,s ? 


Tl‘<“ = -V2g 2 (\Z* + pan OwZ^Zf; - Y'Z 4 \Z< 

2^2 


u* 

i+3,s ? 


-pUiUs _ 

* ■ 3 


-92 ta ne w zt;z% 3i , - , 


(B8) 


into (B7), one recovers the expressions given in [105]. 


To simplify (B7), we expand all mixing matrices in powers of v /Msusy■ At leading order, 


the elements Zf 7 are given by (23), while products of squark mixing matrices simplify as 
follows [165]: 


£ z i. z ik. D ii. = A 

S=1 

£ zizik,D\ h. = A 

S=1 

6 

E yd yd* j-)± _ p± 

Z/ i+ 3 ,s Z/ «+ 3 ,s- Ly ij i q : a 

'Z z i z i D ii. = L i 


■i 9 9 

mz r — rn„ 

'll 1 - 1 


L ± - Rf 

_ CLi CL% 

^di 9 2 


S=1 


S =1 


E yu yu* 7-)± _ p± 

Z/ i+3,s Z/ t-|-3,s- Z - y </ i <j s - r ht; 


S=1 


T z? s z?;d±~ = l± , 


(B9) 


S =1 


where rrZ L and m 2 ;i correspond to the upper and lower diagonal elements of the squark 


(mass) 2 matrices (18), 


A Ui = -m Ui (A™ + /i cot /3) and A di = - m dt (A% + fi tan /3 ), 

are the off-diagonal elements, and the squark propagators in the chiral basis are 

1 


gi — _ 

qi {m x ± m qi ) 2 — m 2 -s + R ’ 


for S = L or R. 


(BIO) 


(Bll) 
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Neglecting terms of 0(u/M| USY ), the final result is 


= -Re | - -g 2 2 t‘a.n6 w {lZf* Z** + ± tan 0 W Z**Z^)A Ui L+.R+ 

- 2 Ag 2 t a . n e w Z£Z£Y u ,R+ + s/2g 2 (±Z£Z£ + itune w Z£Z£)Y^L+ 
+ (L+,R+) ++ (L-.R-) 


91 


= m „.^Re e^ M i 


2 r+p+ 1 + g cot f3) , + + ' 

9^ Ui Ui + 6 |M !| 2 -|^| 2 ^ Ui Ui \ 


+ (L + ,R + ) o ( L" iT), 

"WC'l, = ^Re | tan 6 W Zf* Z**) A di L+ R+ 


v /2 


-V2g 2 (±Z£Z£ - ltan9 w Z%Z£)Y di LX + ^g 2 tan 9 w ZfiZ$Y il B+ 


+ (L + ,R + ) yy ( L~,R~) 
= — fn di ~rRe {e l4>M ' 

+ (L+,R+) O ( L~,R ~), 


l X r+p+ | l (M*+/rtan/ 3 ) + ' 

9 di + 6 |M!| 2 -|p| 2 ^ d ‘ + Kdi \ 


(B12) 


where the squark mixing X q . is defined in (28). 

Higgs exchange —Now consider the Higgs contribution, for which the f-channel ampli¬ 
tude reads 


ax = £ % [*( r 5 p * + r S*n 


L) Ml 


k= 1,2 


-m 


H k J 


s 4 [;(r"$,p R + r^-pj] « 2 


i — 2 ~Re{T^}Re {r^g.}(M 3 Mi)(M 4 -u 2 ) + spin-dependent terms. (B13) 

k= 1,2 


Evidently, the Wilson coefficient due to Higgs exchange is 

c ’fd = 5E^r R *( r S}Re{r«* l }, 


fc=l 


(B14) 


H k 


and as in the case for squark exchange, we substitute the H®xx and couplings [103] 

r"* = — (Z h lk z£ - Z 2 \Z£)(Z£ SW - Z£c w ), 
c w 

Z. 2 k for q = u 


c, = ~\k z * where z * = 


-Z^ k for q = d 


(B15) 
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to obtain the analytic expressions 
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